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Abstract 

We consider a model of F(R) gravity in which exponential and power corrections to Einstein- A 
gravity are included. We show that this model has 4-dimensional Eguchi-Hanson type instanton 
solutions in Euclidean space. We then seek solutions to the five dimensional equations for which 
space-time contains a hypersurface corresponding to the Eguchi-Hanson space. We obtain analytic 
solutions of the F(R) gravitational field equations, and by assuming certain relationships between 
the model parameters and integration constants, find several classes of exact solutions. Finally, 
we investigate the asymptotic behavior of the solutions and compute the second derivative of the 
F(R) function with respect to the Ricci scalar to confirm Dolgov-Kawasaki stability. 
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I. INTRODUCTION 



In recent years, much work has been done in order to bring us more information on extra 
dimensions. In this domain, investigation of black objects, soliton solutions and instan- 
tons are of particular importance, and could play a key role in opening a window to extra 
dimensions. 

In four dimensions, it is well known that some of the static black hole solutions described 
by the Majumdar-Papapetrou solution [I] can be extended to multi-black objects in higher 
dimensions |2] . The four dimensional Riemannian manifolds for gravitational instantons can 
be asymptotically flat, asymptotically locally Euclidean, asymptotically locally flat or com- 
pact without boundary. For example, Hawking's interpretation of the Taub-NUT solution 
[3] is an example of asymptotically locally flat space. The simplest nontrivial example of 
asymptotically locally Euclidean spaces is the metric of Eguchi-Hanson [3] . Asymptotically 
locally Euclidean instantons have been found explicitly by Gibbons and Hawking [5] and 
they are known implicitly through the work of Hitchin [G]. The complex projective space 
CP 2 is an example of compact, anti self-dual instanton solving Einstein's equation with a 
cosmological constant term [7j . 

Some what more recently, new soliton solutions with interesting properties were discov- 
ered a few years ago in 5-dimensional Einstein gravity [9] . These solutions resemble the 
Eguchi-Hanson metrics [I] with AdS/Z p asymptotics. Eguchi-Hanson metrics were origi- 
nally obtained as solutions to the 4-dimensional Euclidean Einstein equations [TT] and have 
the form 

^7^.2 ~,2 

ds 2 = ^ + -((# + cos 6dcp) 2 + d9 2 + sin 2 9dip 2 ) (1) 

1 - 74 

where a is a constant of integration. The Riemann curvature tensor of these solutions is 
self-dual, obeying the relation e a/3 ^ u R^ vpr = R a /3 pT - 

Once we entertain the notion of extra dimensions, we need not impose the asymptotic 
flatness condition, spherical topology of the horizon and other four dimensional restrictions. 
In fact, higher dimensional solutions admit a variety of asymptotic structures and horizon 
topologies that have more interesting properties and richer structure than four-dimensional 
ones. For example, there have been many interesting results in Kaluza-Klein theory [TU] 
with compact extra dimensions as well as in solutions making use of Eguchi-Hanson and 
Taub-NUT spaces [HI [12] with nontrivial topology and structure. 
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From a geometrical point of view, one of the simplest modifications of the gravitational 
interaction to higher order is F(R) gravity, whose action is an arbitrary function of the 
curvature scalar R [TBTfTo] . A key motivation for considering this class of theories has to do 
with addressing a number of cosmological problems. These include the power law early-time 
inflation (TTJ [18] , late-time cosmic accelerated expansion [TEH21] and the singularity problem 
arising in the strong gravity regime [221426] . Furthermore one can find an explanation of the 
hierarchy problem [27J, the four cosmological phases [28] and the rotation curves of spiral 
galaxies [22] within F(R) gravity. 

Following the method of [3D] (which is concerned with a known class of the F(R) theory 
instead of Einstein gravity) the main scope of this work is to present some results on Eguchi- 
Hanson like metrics and investigate their interesting properties. Although many of the 
solutions we obtain have unphysical properties, we find generalizations of the Eguchi-Hanson 
instanton in 4 dimensions, and a new soliton solution in 5 dimensions. 

In contrast with general relativity, F(R) gravity may be considered intrinsically unstable. 
This difference is due to the fact that the field equations in general relativity are of second 
order and therefore their trace gives an algebraic equation for Ricci scalar, but the field 
equations of F(R) gravity are of fourth order and so their trace gives a dynamical equation 
for Ricci scalar. It has been noted that if F(R) ^ and dF/dR = then no stable 
ground state existed [31]. For the solutions we obtain we find that F(R) = dF/dR = 0. 
Furthermore, it has been shown that the effective mass of the dynamical field of Ricci 
scalar is related to the second derivative of the F(R) with respect to R [32] ■ Therefore, in 
order for the dynamical field to be stable its effective mass must be positive, a requirement 
usually referred to as Dolgov-Kawasaki stability criterion. The Dolgov- Kawasaki instability 
criterion, which has found in the metric version of F(R) theories, is sufficiently strong to 
veto some models [33] . 

The outline of our paper is as follows. In section [IT] we present a short review of field 
equation of d-dimensional F(R) gravity. This field equation is then solved for 4-dimensional 



Eguchi-Hanson space in Sec. Ill We then generalize it to 5-dimensional Eguchi-Hanson like 
space-time and obtain a new soliton solution. Conclusions are drawn in the last section, and 
we present more solutions using the same approach in the appendix. 
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II. BASIC FIELD EQUATIONS 



We start from the following action of pure general F(R) gravity 




Variation of Eq. ^ with respect to metric g^ u , leads to the field equation of F(R) gravity 





where R^ u is the Ricci tensor and F R = dF(R)/dR. We consider a recently proposed [5H 

F(R) model 



in which exponential and power corrections to Einstein-A gravity are included. This choice 
of F(R) gravity has some interesting properties such as providing charged solutions from 
pure F(R) gravity [30J. 

Viable modifications to Einstein gravity must pass all sorts of empirical tests, from the 
large scale structure of the universe to galaxy and cluster dynamics to solar system tests. 
One of the outstanding questions in F(R) gravity is whether it is consistent with solar 
system tests or not. When the correction term to Einstein gravity is of the exponential form 
[T4"| 134"] . it can be shown that there is no conflict with solar system tests and that the high 
curvature condition is satisfied [35]. Furthermore, the solutions of this model are virtually 
indistinguishable from those in general relativity, up to a change in Newton's constant [35J. 

In addition, one can choose R 2 with exponential corrections to Einstein gravity [54"l I5B] to 
explain inflation. By adjusting the free parameters, this model can satisfy the high curvature 
condition, early universe inflation, stability and the late time acceleration [35]. Also, in order 
to resolve the singularity problem arising in the strong gravity regime, Kobayashi and Maeda 
[26J have considered a higher curvature correction proportional to R n where n > 1. 

Motivated by the above, we seek instanton solutions to F(R) gravity using Q for 4- 
dimensional Eguchi-Hanson space and soliton solutions for Eguchi-Hanson like (4 + 1) - 
dimensional space-time. 



F(R) = R-2A-X exp(-£#) + kR 



•n 



(4) 
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III. GENERALIZED EGUCHI-H ANSON INSTANTONS 



We consider first the 4-dimensional Eguchi-Hanson type metric ansatz 

ds 2 = ^ + r 2 (al + a}) + r^r)** (5) 
where the differential one forms Cj can expressed in terms of Euler angles 9, <fi and i/j 

&x = - (sin ipd9 — sin 9 cos ipdip) , 
o y = - (— cos ipd9 — sin 9 sin -^g^) , 

(dip + cos 0d<^) . (6) 
For F(R) = R, it is well known that the solution to the Euclidean Einstein equations is [I] 

9<r) = i - £ (r) 

where a is constant. For a = the constant-r surfaces of the metric ^ have the geometry 
S 3 , whereas for a/0 the geometry must 5 ,3 /Z 2 to avoid singularities. Note that the metric 
a x + a y yields the geometry of a 2-sphere of radius 1/2. 

We wish to obtain general solutions using the ansatz ^ for the mentioned model of F(R) 
theory, and investigate their geometrical properties. 

Since the metric is purely Euclidean with no time-like direction, any solutions we obtain 
cannot be interpreted as black holes. If the function g(r) vanishes at r = r then the metric 
will have a singularity unless ip is appropriately periodically identified with period 

a , 167T 

M = Ti — \ 8 

where g' = dg/dr. 

Using Eqs. ^ and ([3]), we find the solution 

2 2 

A = -fie#, (10) 

for the metric function, where II = n/tx n_1 + 1, f2 = 2A — nx n ~ X-, an d X> Q an d a are 
integration constants and £, /3, k and n are free parameters. In order to study the geometrical 




FIG. 1: Eq. g(r) versus r for \ = — 1, q = 1, and a = — 1 (solid line) and a = 1 (dashed line). 



structure of this solution, we first look for the essential singularity(ies). The Ricci scalar 
and the Kretschmann scalar are 

R = x, (ii) 

and 



1 L fj,upa- 1 L 



384a 2 384ag 2 128g 4 x 2 



(12) 



which means that the Kretschmann scalar (12) diverges at r = 0, is finite for r / and 
approaches y as r — > oo. 

For x < and positive a this solution has a singularity at r = 0. However for \ < and 
negative a the function g(r) vanishes at a finite value of r we denote by r . The space is 
geodesically complete provided the coordinate ip nas period 

^o 2 



where 



3|xl r o + 8<? 2 + 16rQ 



24 



+ + 1 2r o 







Figure [I] illustrates that for negative a, the function g(r) vanishes at a certain value of r; 
this solution is like a "bubble" in Euclidean space or a kind of soliton. 

For x > and positive a the function g{r) has a real root at Tq, becoming negative for 
larger values of r. This space has a naked singularity at r = 0. However if a is also negative 




FIG. 2: Eq. g(r) versus r for % = 1, g = 1, and a = — 1 (solid line) and a = 1 (dashed line). 

then the metric function g(r) has two roots, r< and r>, which delineate the range of r (see 
Fig. [2j. The solution will be free of singularities provided the period of ip at r< is the same 
as that at r>. We find that there are no real values of the parameters for which this situation 
holds, and so conclude that all x > solutions are singular. 

This solution is problematic insofar as the circle described by ip is growing exponentially 
relative to the 2-sphere as r — > 00. The reason is that g(r) is growing like r 2 for large r, so 
the ^-circle is growing much faster than the 2-sphere for large r. This will not happen if 
X = 0, which implies £ = and A = — 2A. 

We therefore conclude that x — and a < yields the only physically reasonable class 
of instanton solutions, whose metric is 



ds 2 



dr 7 



1 4- 2- 



^+r 2 K 2 + aJ)+r 2 1 + 



Si - °L - 

r 2 r 4 J ^z' 



(13) 



where we have set a = —a 4 . We find that the metric function vanishes at r 

V_2( v / g 4 + 4a 4 _ g 2) and that 



Aif) = 4tt 1 



(14) 



generalizing the Eguchi-Hanson metric ^ to -F(-R) gravity. Note that the class of metrics 
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FIG. 3: Eq. (15): Frr versus K for \ = ~~ 1> A = 1 and n = 3. 



(13) cannot be obtained by Wick rotation of some coordinate in a corresponding Lorentzian- 



signature space-time. 

Although we have not included a time-like coordinate yet, the F(R) model can be checked 
for stability. In order to check the stability condition (Dolgov-Kawasaki stability), we obtain 
the second derivative of the F(R) function with respect to Ricci scalar 

n 2 



Frr = ^r + n(n- l)«x 



n-2 



(15) 



which shows that this model is stable provided the free parameters of the model are chosen 
suitably (see Fig. [3] for more details). We should note that for x — 0, Frr reduces to 

1 



RR 



2 A' 



(16) 



which confirms that this solution is stable for A > 0. 



IV. 5-DIMENSIONAL SOLITON SOLUTIONS 



Here, we include the time coordinate and begin with the following 5-dimensional ansatz 
for the metric 



ds 2 = -h(r)dt 2 + 



dr 2 
h(r)g(r] 



+ r 2 K 2 + trj) + r 2 g(r)o\ 



(17) 
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In this metric, the hypersurface t =constant corresponds to Eguchi-Hanson space when we 
fix h(r) = 1 and set g(r) as in Q. For g(r) = 1, Eq. ( fl7"| ) reduces to 

*-7/v>2 ^2 ~,2 

rfs 2 = -h{r)dt 2 + — + — (d# 2 + # 2 + # 2 ) + — cos(0)##, (18) 

which is Schwarzschild-like space-time. When we choose h(r) = 1, one can find that g{r) 
is the same as Eq. Q with the same analysis. In this section, we are looking for exact 
solutions of the field equations ^ for the Eguchi-Hanson like metric ansatz (17) with F(R) 
given by Eq. Q. 



We should note that for any possible solution for the metric ansatz (17), one may en- 
counter several cases of asymptotic structures and singularities of the metric: 

• h(r) flips sign at one or more values of r and g(r) does not change sign for sufficiently 
large r. These correspond to the usual kinds of horizons: black hole and cosmological. 
If h(r) — > 1 at large r then there will be no cosmological horizon. There will be a 
black hole if h(r) flips sign at small r before g(r) does. 

• g(r) flips sign at one or more values of r and h(r) does not change sign. At small finite 
r > this will be the edge of the soliton (or bubble), and if g(r) — > 1 at large r then 
the space-time will be a soliton. If g(r) flips sign again for large r then the space will 
be compact. 

• g(r) flips sign at one or more values of r and h(r) changes sign at other values of r. If 
g(r + ) = and h(r c ) = then this is a soliton surrounded by a cosmological horizon. 
If h(r) has multiple roots, all larger than r + , then the soliton is in some kind of black 
hole that has multiple horizons; any roots in the range r < r + don't matter because 
the soliton covers them up. 

• g(r) and h(r) both flip sign at the same value of r. In this case the coordinate 
ip becomes time-like and the space-time has closed time-like curves (CTCs) beyond 
this value of r (either larger or smaller, depending on where g(r) and h(r) are both 
positive). 

The only way to get a sensible AdS, or dS, flat structure asymptotically is if g(r) — > 1 
and h(r) — > ±r 2 or h(r) — > 1 for large r - anything else gives some other geometry. For 
example if h(r) — > ±r p (and g(r) — > 1), then the geometry will be asymptotic to a Lifshitz 
geometry. 



1. First solution set: 



We find a more interesting class of space-times if we set A = e€x [3x 10 ^ x lf n ^ x - ■ Solving 
the field equations ^ and (17) we find 

t 

.-A 



9{r) = 1-^, (10) 



h(r) = 1 — — r 2 . 
v ; 20 

After some algebraic manipulation, we find that the Kretschmann and the Ricci scalars 

are 

o „ _ 384a« _ 96^ 9aV , £ (9n , 

^ ' r 12 5r 10 + 50r8 + 10' 1 J 

R = X, (21) 

indicating a curvature singularity at r = 0, and a space-time of constant curvature as 



r — > oo. The singularity in (20) is not accessible from the space-time, which is free of 
singularities provided the period of ip is appropriately chosen [8]. This solution is the same 
as the Eguchi-Hanson soliton [8J, |9], and will therefore share all of its properties. 

To investigate Dolgov-Kawasaki stability, we compute derivative of the F(R) function 



RR 



[«(2n - 5)x n - 3x + 10A] £ 



2 



+ n(n-l) KX n ~ 2 , (22) 



2^ + 5 

and see that for suitable values of free parameters this model is stable (see Fig. [4]). 
2. Second solution set: 

Here, we produce another set of solutions with some limitation on the free parameters 
of the model. We will see that for a certain choice of parameters we obtain a new soliton 
solution that generalizes the Eguchi-Hanson soliton. 

Choosing the following parameters 

n = 0, 

A = e«+x-2A ( K + x - 2A) , 

( ^ (^2A)< (23) 
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FIG. 4: Eq. (22): F RR versus re for x = -1, A = 1, n = 3 and £ = 1. 



we find a new solution set of the field equation ^ with metric function (17) 



a ( r \ = i _ 5Mx _ xif 

y ^ > 96 24 r 3vP=M' / 2 4) 

/i(r) = l M 



r 2 ; 



where M, x an d & are integration constants, and we have redefined k — 2A as an effective 
cosmological constant h- e ff- 

Unless x = t ne metric will not have the desired asymptotic behavior, and so we make 
this choice, yielding an asymptotically (locally) flat metric. Note that for positive M this 
solution is real for all r > M; at r = M it is singular. However this singularity is excised 
from the space-time for all b > 0, since the function g{r) has a single root at r > M, and 
this root is larger than the positive root of h(r), as illustrated in figure [5j 

The metric is 



= - (i - g) + f_ + n W + + ^ (i - 



r 3 \/r 2 — M 



(1 - f ) (l 

(25) 

and is an interesting new generalization of the EH soliton. The solution is nonsingular and 
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FIG. 5: Eq. (24): g(r) (solid line) and h{r) (dashed line) versus r for x = 0, M = 1 and b = 0.5. 



the quantity M can be of either sign (see figure [6] for the behavior of the metric functions 
for M < 0), with the size of the soliton an increasing function of M. 

The soliton radius r s is given by solving a quartic equation, which yields 



r s ± 



V\M\ 



sgn(.U) ( - + — ^9 + 6yV3 _ 288<T 8 y-V3 



+ 



12 



/ 3 _yi/3 + 48(T 8y-i/3 



9 



-.1/2 



V9 + 6yV3 _ 288a 8 r-V3 



(26) 



where = 6 4 /M 2 and F = 12aV768(x 8 + 81. Each of r s ± are increasing functions of cr, 
growing linearly with a for a sufficiently large. 

To ensure regularity of the soliton, the periodicity of ip must be chosen so that 



Aijj 



(27) 



(4r2 - 3M) p 

where p is an integer, the latter equality following from the elimination of string singularities 
at the poles. Note that if M — then p — 2, consistent with the Eguchi-Hanson instanton. 

The above relation constrains the value of a for any integer p. It is straightforward 
to show that solutions exist for all integer values of p > 1. If M > then the solution 



12 



FIG. 6: Eq. (24): g(r) (solid line) and h{r) (dashed line) versus r for x = 0, M = —1 and 
b = -0.5. 

a = 0.720202 for p = 2, decreasing as p increases, whereas if M < then the solution 
cr = 1.090507733 for p = 2, decreasing toward unity as p increases. 

The ADM mass of this solution can be calculated using, for example, the background 
subtraction method, and is easily shown to be M. = 3iiM/8p. 

3. Third solution set: 



Here we present second solution of Eqs. ^ and (17) with R = \ as 



9{r) 
h(r) 

where we should adjust 



-3x a 4 



+ 



10C r 4 rVCr 2 + 6 



+ 



4(10C + 3 X ) 
5C 2 r 2 



, C 2 
o 



(28) 



K = 0, 

A = ( x _ 2A) 
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which confirm that considering Eq. (|28j), leads to vanishing the R n correction. For this 
solution to have physically reasonable asymptotic properties, we must set x — — 10C/3, 
yielding 



a 4 

g(r) = 1 



b 



r 4 r4v / r 2/£2 + 1 ' 
2 

h(r) = 1 + ^, (30) 



where to have real solutions, we have set C = 6/£ 2 > 0, and we have rescaled b — > b for 
convenience. If b > then the metric function g(r) has one real root at r = r s > 0, given 
by the square root of the real solution to the equation 



which is a quartic equation in y/r~ s . For large r this space-time is asymptotic to the Eguchi- 
Hanson soliton [8] . The curvature diverges at r = (and hence not at a point located within 
the space-time). For large values of r, one obtains 

lim R a ^R a ^ = ^. (31) 

Regularity in the (r, if)) section and also elimination of string singularities at the north 
and south poles implies that 

5r 4 + Ar 2 s l 2 - a 4 ~ p ^ ' 

Here, we analyze the Dolgov-Kawasaki stability. It is easy to show that 

I 2 



RR 



20 + 2A/ 2 



(33) 



and therefore we conclude that this solution is stable for A > — 

It is straightforward to show using the methods of Ref. [3Zj that the mass of this solution 
is M. = —7ra A /8£ 2 p, the same value as for the Eguchi-Hanson soliton |8J, apart from a 
constant factor due to the Casimir energy. 
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V. CONCLUSIONS 



In this work, we have considered a kind of well-known F(R) gravity in Eguchi-Hanson 
space and Eguchi-Hanson like space-time. We have shown that, in 4-dimensional Eguchi- 
Hanson type Euclidean metric, the solutions can be interpreted as stable instantons. 

Upon including the time coordinate we found two distinct generalized 5-dimensional 
Eguchi-Hanson space-times in F(R) gravity, one that is asymptotically flat and another 
that is asymptotically AdS. These are obtained by imposing additional constraints on the 
F(R) model parameters. We have investigated their asymptotic behavior, computed their 
masses, and obtained constraints on the parameters to ensure regularity of the metric. Our 
investigation of Dolgov-Kawasaki stability of these solutions indicates that their objects are 
stable. 



It is possible to obtain more solution sets for Eguchi-Hanson like (17) space-times using 
our approach. However an inspection of their basic properties indicates that they are not 
physical (see appendix for more details). 

It would be worthwhile to investigate the thermodynamic as well as dynamical stability 
of solutions we have found. Other interesting questions, such as whether or not the sub-class 
of theories yielding the EH-like metrics we have found obey Birkhoff's theorem |38j, remain 
interesting problems for future consideration. 
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APPENDIX 

A2: Fourth and fifth solution sets: 

We produce here two solution sets with the same limitation on the free parameters. 
Considering the model parameters in the third solution set, one can find two additional 
solution sets of the field equation ^ with metric function (17) as follows 




. =l-Cr, 

IV '■ \ . 2(41C 2 -12 X ) _ 6(205C 2 -8x) {2 _ ( 34 ) 

41C r 533C 2 2665C* 3 \r Cr 2 ) ' 

a (r) - ^ r + - 8 ( 82c2 ~ 5x ) (1 + -*-) 

y . I y\' ) 41C*' ^ 82C 2 451C 3 \r ^ Cr 2 ! i 

h(r) = 1 — Cr, 

where x and C are integration constant. Since n — 0, and as we noted for third solution 
set, one can redefine k — 2 A as an effective cosmological constant. 

Calculating the Kretschmann scalar leads to the following asymptotics 

lim^R^R^ — »• oo, 

lim R a ^ v R a ^ v = (36) 

r — >oo loot 

Considering solution III with negative C, one can find that h(r) flips sign at large r and 
g{r) doesn't, and so we have a cosmological horizon. Respectively adjusting x to ^jff- and 
52 4i9 72 * n solutions III and IV, we find that these space-times have a causal horizon since the 
functions g(r) and h(r) both have the same root in r = \jC. 

One can note that these kinds of solutions both have the unphysical asymptotic structure. 

Here, we compute second derivative of the F(R) function for the mentioned three sets 

-1 



1 RR 



K + x-lk 1 

which confirm that this solution is stable for A > | (k + %). 



(37) 



A3: Sixth and seventh solution sets: 

Looking at Eqs. (|3]) and (17), and considering Eq. (10), one can obviously obtain two 



solution sets with constant curvature scalar (R = x) as follows 

9(r) = 1 - £ r 2 , 



VI : < b , , (38) 

' - 3 X~C I 54( X -4C) v 1 

b V I 11C ^ 55C 2 r 2 ' 



VII : { 9{r) 22 ^ 2 5 ^ 4 ' (39) 
h(r) = 1 - Cr 4 , 
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where k and n are free. 



For the sixth solution set (Eq. (38)), calculations show that the Kretschmann scalar 
diverges at r = 0. It is more interesting to note that notwithstanding previous solutions, it 
depends on two parameters as r — > oo 

Hence there is a curvature singularity at r = 0. This singularity will be hidden behind an 
event horizon provided either 

4C > X > j 

provided C and \ are both positive. However the space-time ends at the coordinate value 
r = a/6/C, and so does not describe a black hole with familiar asymptotic properties. 
For C and x both negative there will be an event horizon cloaking the singularity provided 

\X\ >4|C| 

In this case the space-time geometrically consists of a squashed 3-sphere whose squashing 
parameter grows exponentially with r. 



The seventh solution set (Eq. (39)) has a curvature singularity at r = 0; for large r, we 
obtain 

lim R a( s^R a ^ = (41) 

r — >oo 121 

Note that the circle described by ip is growing exponentially relative to the 2-sphere as 
r — > oo and so these solutions both have the same types of problems as the instanton 
solutions. In order to check the Dolgov-Kawasaki stability, we should obtain Frr. It is easy 



to show that F RR is the same as Eq. (15). 
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